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ABSTRACT 

We present the new constraints on the cosmic optical background (COB) obtained from an analysis 
of the Pioneer 10/11 Imaging Photopolarimeter (IPP) data. After careful examination of data quality, 
the usable measurements free from the zodiacal light are integrated into sky maps at the blue (^0.44 
/im) and red (^0.64 /mi) bands. Accurate starlight subtraction is achieved by referring to all-sky 
star catalogs and a Galactic stellar population synthesis model down to 32.0 mag. We find that the 
residual light is separated into two components: one component shows a clear correlation with thermal 
100 fim brightness, while another betrays a constant level in the lowest 100 /im brightness region. 
Presence of the second component is significant after all the uncertainties and possible residual light 
in the Galaxy are taken into account, thus it most likely has the extragalactic origin (i.e., the COB). 
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The derived COB brightness is (1.8 ± 0.9) x 10~ 9 and (1.2 ± 0.9) x 10 

the blue and red band, respectively, or 7.9 ± 4.0 and 7.7 ± 5.8 nW m~ 2 sr _i . Based on a comparison 
with the integrated brightness of galaxies, we conclude that the bulk of the COB is comprised of 
normal galaxies which have already been resolved by the current deepest observations. There seems 
to be little room for contributions of other populations including "first stars" at these wavelengths. 
On the other hand, the first component of the IPP residual light represents the diffuse Galactic light 
(DGL) — scattered starlight by the interstellar dust. We derive the mean DGL-to-100 /im brightness 
ratios of 2.1 x 10 -3 and 4.6 x 10~ 3 at the two bands, which are roughly consistent with the previous 
observations toward denser dust regions. Extended red emission in the diffuse interstellar medium is 
also confirmed. 

Subject headings: cosmic background radiation — cosmology: observations — dark ages, reionization, 
first stars — dust, extinction — galaxies: evolution — infrared: ISM 
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1. INTRODUCTION 

The cosmic optical background (COB) is the optical 
component of the extragalactic background light, which 
is the integrated radiation from all light sources out- 
side the Galaxy. Dominant contribution to the COB 
comes from stellar nucleosynthesis in galaxies at redshifts 
z < 10, while other mechanisms such as mass accre- 
tion to super massive black holes in active galactic nuclei 
(AGNs), gravitational collapse of stars, and particle de- 
cay can contribute to the COB brightness. As a fossil 
record of light production activity in the Universe, the 
COB conveys information on the cosmic star formation 
history including birth and death of Population III stars. 
Comparison of the measured COB brightness and the in- 
tegrated brightness of resolved galaxies has fundamental 
importance, since it provides a direct clue as to the frac- 
tion of celestial visible light whose origin is known to us. 

However, robust detection of the COB has long been 
hampered by the extremely bright foreground emissions. 
While expected brightness of the COB is around 1 bgu 
= 1 x 10~ 9 erg s" 1 cm -2 sr _1 A -1 , the terrestrial air- 
glow and the zodiacal light (ZL) are a few orders of mag- 
nitude brighter tha n this level at optical wavelengths 
(|Leinert et alj|j"998l ). Besides these brightest two com- 
ponents, the diffuse Galactic light (DGL), which refers 
to scattered starlight by the interstellar dust, is another 
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but much fainter component of t he diffuse light of the 
night sky. iDube et al.l (|1977lll979f ) attempt a COB mea- 
surement by subtracting all the foreground components 
from the night-sky brightness observed at Kitt Peak Na- 
tional Observatory. The terrestrial airglow is assumed to 
follow a function of the zenith angle with temporal varia- 
tion, while the ZL brightness is estimated from equivalent 
width of an imprinted solar Fraunhofer line on the ZL 
spectrum, making use of the fact that the solar spectrum 
is well reproduced in the ZL due to weak wavelength 
depe ndence in scatt ering strength of the interplanetary 
dust. iMattilal (|1976l ) devised the method to utilize Galac- 
tic dark clouds as shielding walls of the COB. By com- 
paring the brightness measured at two lines of sight, one 
toward a dark cloud and one toward a nearby "trans- 
parent" sky, the background light coming from behind 
the dark cloud could be extracted. However, these and 
other similar attempts were not very successful due to 
inadequate knowledge about the foreground radiations, 
especially spatial and temporal variations of their bright- 
ness. 

One of the fundamental solutions to the problem of 
foreground subtraction is to carry out measuremen ts at 
outside the relevant emission regions. iToIierl CL983) uti- 
lizes the measurements of the night-sky brightness by 
a photopolarimeter on board the Pioneer 10 spacecraft. 
He analyzes the data obtained at the heliocentric dis- 
tances larg er than 3.3 AU, wh ere the ZL is found to be 
negligible (Ha nner et al.lli"974l ). Unfortunately, though, 
Galactic stars cannot be resolved due to poor spatial res- 
olution of the instrument, which dominate the measured 
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sky brightness and prevent robust detect ion of the COB. 

The latest attempt is made by iBernstein et alj 
( 2002a. b.c), who carried out a coordinated observation 
program using the Hubble Space Telescope (HST) and a 
ground-based telescope at Las Campanas Observatory. 
The terrestrial airglow is avoided by using space obser- 
vations, w hile the ZL b r ightness is e stimated by a similar 
method to lDube et al.1 (|1977l [l979h . i.e., from measured 
equivalent width of a solar Fraunhofer line reproduced in 
the ZL. The DGL brightness is estimated using a simple 
scattering model and thermal emission brightness of the 
interstellar dust. Ho wever, a number of problems in their 
analysis are raised bv lMattilal ([20031 ). e.g., inappropriate 
assumptions in deriving the ZL and DGL brightness. A 
se ries of corrections to th e ori ginal pape r s are presented 
in IBernstein et all ([20051 ) and IBernstein! ([20071 ). and the 
derived COB brightness is now accompanied by signifi- 
cant errors; the latest values are (6 ± 4), (10 ± 5), and 
(7 ± 4) bgu at the 0.30, 0.55, and 0.80 /im band, respec- 
tively. 

Here we present the new constraints on the COB 
brightness obtained from an original analysis of the 
night-sky data taken by the Imaging Photopolarimeters 
(IPPs) on board the Pioneer 10 and 11 spacecrafts. The 
present work giv es a significa nt improvement of the pre- 
vious analysis by iTollern) 19831 ) by taking advantage of the 
additional higher-quality (Pioneer 11) data, accurate all- 
sky star catalogs, and a realistic DGL model combined 
with the far-infrared (IR) space observations. This paper 
is organized as follows. In Section 2 we describe the ob- 
servations and reductions of the data used in this work, 
including data assessment and map creation processes. 
Subtraction of the Galactic light, namely, the starlight 
and the DGL, are presented in Section 3. Then we show 
the results in Section 4, followed by discussion on the de- 
rived correlation between the DGL and the far-IR emis- 
sion and on the resolved fractions of the COB in Section 
5. Finally a summary appears in Section 6. Hereafter 
the Pioneer 10 and 11 spacecrafts are sometimes abbre- 
viated as "P10" and "Pll" . All the magnitudes are given 
in the Vega-based system. 

2. OBSERVATIONS AND REDUCTIONS 

Here we describe the observations and reductions used 
in this work. While a similar analysis is carried out by 
IGordon et al.l (|1998l) for the investigation of the Galac- 
tic diffuse emission, we clarify the two major differences 
between their and the present analysis as follows. (1) 
Stellar contribution is removed from the each instru- 
ment field of view (FOV) before creating the final, high- 
resolution sky maps in this work, rather than first creat- 
ing star-contaminated maps and then evaluating the stel- 
lar contribution in them. Considering the large instru- 
ment FOV (low angular resolution), the former method 
is more suitable than the latter to derive the spatial dis- 
tribution of the diffuse emission as we will see below. (2) 
The data quality is evaluated at a few stages in the mid- 
dle of the data processing in this work, which allows us 
the accurate understanding of the uncertainties in the 
final results. 

2.1. IPP Observations 

The Pioneer 10 and its sister ship Pioneer 11 are 
NASA's space missions, which were launched on 1972 



March 2 and 1973 April 5, respectively. They are the 
first spacecrafts to travel through the Asteroid belt and 
explore the outer solar system. One of the scientific in- 
struments on board the spacecrafts is the IPP, which 
was designed to make two-color (blue and red) maps 
of the brightness and polarizati on over the sky dur - 
ing cruise phase of the missions (IPellicori et al.l 119731) . 
The IPP data taken at the heliocentric distances R — 
1-5 AU were used to study distribution of the inter- 
planetary dust responsible for the ZL an d diffuse light 
components outside the ZL clouds (e.g. . lHanner et al.l 
\W74 . IWeinberg et al.l Il97l iToller et al.1 119871 ). Below 
we give a description of the IPP measurements and the 
primary data p rocessing. For more det ails, s ee, e.g., 
IWeinberg eFax] (| 19741 ) and IGordon et al.1 (119981 ). 

The IPP is equipped with a 2.5 cm Maksutov-type 
telescope, together with a Wollaston prism, multi- 
layer filters, and two dual-channel Bendix channeltrons 
(IPellicori et al.|[l97l . It measures simultaneously two 
orthogonal polarization components of brightness in two 
(blue at the wavelengths 0.39 - 0.50 /im and red at 0.60 
- 0.72 pm; "5i PP " and "i?i PP " hereafter) broad-band fil- 
ters. Effectiv e wavelengths of the two bands are 0.44 /jm 
and 0.64 /mi (jGordon et al.l[l998h . The two polarization 
components are defined with "S" vector parallel to the 
instrument rotation axis and " P" vector which is perpen- 
dicular to it. Thus one IPP exposure generates four sur- 
face brightness measurements, S (S-vector component) 
and S p (P-vector component) at each of the i?i PP and 
i?i PP band. The four measurements are carried out by 
the different photocathodes placed at the end of the dif- 
ferent optical paths. The data used in the present work 
was obtained with the " ZL mode" , with instantaneous 
FOV of 2°. 29 x 2°. 29. Sky brightness is integrated for 
1/64 (one "sector") of a 12.5 s spacecraft spin period, re- 
sulting in effective FOV of 2°. 29 x (2°. 29 + 5°. 625 sini) 
where L is the angle between the instrument pointing 
and the spacecraft spin axis (called " look angle" ) . With 
a change of look angle, the effective FOV ranges from 7.5 
deg 2 (L = 170°) to 18 deg 2 (L = 90°). Since the single- 
pixel detector is used in each of the four channels, the 
FOV of the detector size also gives the angular resolution 
of the instrument. A data cycle of the IPP measurements 
consists of 10 rolls; eight for sky measurements, one for 
photometric calibration using a self-luminous phosphor 
source, and one for monitoring offset and dark-current 
levels. Sky maps are created by changing look angle L 
and repeating the data cycles. 

Primary processing of the obtained data has been car- 
ried out by the IPP team. Initially absolute photomet- 
ric calibration is achieved using the pre-flight measure- 
ments of a 14 C-activated self-luminous phosphor source. 
Instrument pointings (coordinates of FOV centers) are 
computed from the spacecraft attitude parameters. The 
brightness contribution of bright stars is subtra cted re- 
ferrin g to stars in the Yale Bright Star Catalog (Hofficit 
H96l and V < 8 stars in the U.S. Naval Observ atory 
(USNO) Photoelectric Catalog ([Blanco et al.lll968| ). Fi- 
nally, observed star crossings and in-flight measurements 
of the internal phosphor source are used to correct for 
temporal variation of the instrument sensitivity and to 
revise telescope pointings. The final position accuracy 
of FOVs is 0°.15 - 0°.40. The processed data are stored 
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in th e Background Sky Tape (Weinber g fe Schuerman I 
[l98Th . which is now available from NASA's National 
Space Science Data Center. 

Using the Pioneer 10 IPP data collected at the var- 
ious heliocentric distances R, lHanner et al.1 (|1974l ) find 
that the ZL brightness at R = 2.41 AU is less than 10 % 
of that observed at R = 1 AU, and that the ZL bright- 
ness is below the detectable level of the instrument when 
the spacecraft is at or beyond R = 3.26 AU. Hence, the 
data obtained at larger distances are most suitable for 
an analyses of diffuse radiations outside the ZL clouds. 
In Table [1] we summarize the IPP measurements carried 
out at R > 3.26 AU. 

2.2. Data Assessment 

First, we look through all the measurements in Ta- 
ble [1] in order to identify and remove the poor-quality 
data. Recall that one IPP exposure generates four 
measurements, i.e., S p and S s at each of i?ipp and 
i?ipp. We remove the exposures with a negative flux 
record in at least one of the four measurements, which 
are most likely the corrupted data due to interrupted 
data communication between the spacecrafts and the 
ground station. Also removed are the exposures taken 
when the angle between the instrument pointing direc- 
tion and the sun ("elongation") is less than 70° for the 
Pioneer 10 and 45° for the Pioneer 11. These expo- 
sures could be subje ct to the scattered sunlight inc ident 
on the instruments (| Weinberg fc Schuerman I II 98 lh . We 
then calculate median and standard deviation of the four 



and 



flux ratios, 5g Ipp /S| Ip 

S% I S\ , of 64 exposures constituting one data roll, 
and remove all the exposures in the data roll showing the 
outlying values. In Figure [1] we show the distributions 
5 b,pp> SjLpi and s \ vv brightness for the 



of the S% . 

iJlPP > 



P10 and Pll measurements after the above removal of 
the poor-quality data. Overall they show the similar dis- 
tributions to each other, while that of the P10 Sjj is 
significantly different from others. It is due to the excep- 
tionally large noise in this channel, therefore we decide 
to discard all the P10 S% data at this stage. 

Kipp o 

The dominant contributor to measured IPP brightness 
is Galactic stars. Since integrated starlight over an FOV 
of ~10 deg 2 is expected to polarize very little on aver- 
age, the two orthogonal polarization components S p and 
S have essentially the equal brightness. Therefore the 
comparison between the two, which were measured by 
the different detector units at the end of the different 
optical paths, can be used to estimate accuracy of the 
IPP measurements. We plot in Figure [2] distributions of 
the brightness deviation 



D = 



S p -S s _S P - S/2 _ S s - 5/2 



S p + S s 



S/2 



S/2 



where S — S p + S s is the total sky brightness at -Bipp 
or i?ipp, in the brightness range S = 30 - 90 bgu where 
this work focuses on. Narrowing the range to S = 30 - 
50, 50 - 70, or 70 - 90 bgu does not change the following 
estimates significantly. Figure [5] shows that the distribu- 
tions of D are close to the Gaussian distribution, there- 
fore mean and standard deviation of the best-fit Gaussian 
functions give an appropriate estimate of systematic and 
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Sg ]pp (panel (a)) and 
(panel (b)) brightness for the P10 (black lines) and 



Pll (gray lines) measurements. Those of the Sg 
shifted upward by +200 for visibility. 
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Fig. 2. — Distributions of the brightness deviation D = (S F — 
S S )/(S P + S S ) in the P10 Bi PP (panel (a)), Pll B WP (panel (b)), 
and Pll .Ripp (panel (d)) measurements (histograms). The dotted 
lines are drawn at D = 0. The superposed smooth curves are the 
Gaussian distributions fitted to the data. 

random uncertainties inherent in the IPP measurements. 
Thus, obtained values of the fractional systematic and 
random errors are 3% and 6% at -Bipp and 3% and 4% 
at i?ipp, respectively, relative to S p or S s . These val- 
ues are divided by 2-^/53 to give the corresponding errors 
relative to the total sky brightness S — S p + S . The 
derived estimates are summarized in Table [5] (first row) 
along with other estimates derived below. 

As already introduced, the total sky brightness is the 
sum of the two polarization components: S — S p + S s . 
We substitute the P10 S^ fluxes for the discarded 

i-r.iPP 

3 The factor 2 accounts for the change to the relative values to 
S, from relative to S p or S s , while the factor \f2 comes from the 
reduction of errors as the independent measurements S p and 5 s 
are summed. 
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TABLE 1 

Observation Journal of the Pioneer 10/11 IPPs at Outside the ZL Clouds 







H a 


Spin Axis b 








Number of 


Usable 


Spacecraft 


Day/Year 


(AU) 


(R.A., Decl.) 


Look Angle 


Exposures 


Fraction 


Pioneer 10 


354/1972 


3.26 


(108° .83, 23° .62) 


30°. 75 - 


170° 


.11 


5,696 


60.3 % 


Pioneer 10 


093/1973 


3.92 


(138° .50, 17° .54) 


128° .62 - 


- 169° 


.98 


1,344 


75.4 % 


Pioneer 10 


149/1973 


4.22 


(141°. 99, 16° .68) 


28° .63 - 


170° 


.66 


5,504 


61.4 % 


Pioneer 10 


216/1973 


4.54 


(135° .03, 19°. 70) 
(130° .60, 19°. 71) 
(127°. 40, 20° .25) 


91° .66 - 


170° 


.07 


2,816 


50.2 % 


Pioneer 10 


237/1973 


4.64 


29° .00 - 


169° 


.40 


5,312 


56.6 % 


Pioneer 10 


279/1973 


4.81 


28° .92 - 


169° 


.22 


5,248 


57.4 % 


Pioneer 10 


021/1974 


5.08 


(143° .97, 15° .21) 


30° .64 - 


169° 


.51 


5,504 


58.4 % 


Pioneer 11 


057/1974 


3.50 


(155°. 13, 11° .95) 


29° .65 - 


169° 


.00 


4,544 


72.6 % 


Pioneer 10 


068/1974 


5.15 


(155° .88, 10°. 23) 


28° .92 - 


169° 


.87 


5,376 


67.2 % 


Pioneer 11 


106/1974 


3.81 


(166° .08, 7° .67) 


29° .63 - 


170° 


.64 


4,672 


73.5 % 


Pioneer 11 


148/1974 


4.06 


(173°. 70, 4° .63) 


29°. 56 - 


170° 


.65 


4,672 


77.1 % 


Pioneer 11 


178/1974 


4.22 


(175°. 11, 4°. 12) 


29°. 54 - 


170° 


.55 


4,608 


72.9 % 


Pioneer 11 


236/1974 


4.51 


(171°. 01, 5°. 94) 


29°. 61 - 


169° 


.93 


4,544 


80.8 % 


Pioneer 11 


267/1974 


4.66 


(165°. 58, 8°. 00) 


31°. 45 - 


168° 


.39 


4,416 


70.4 % 



Note. — a Heliocentric distance of the spacecraft. b The end of the spin axis pointing toward the Earth. 



TABLE 2 

Error Estimates of the IPP Measurements 



Subset Criteria 


Bipp 
Systematic 


Random 


i?ipp 
Systematic 


Random 


Polarization Component 


0.01 


0.02 


0.01 


0.01 


Spacecraft 


0.01 


0.02 


0.001 


0.02 


Heliocentric Distance 


0.008 


0.02 


0.0001 


0.02 


Elongation 


0.003 


0.02 


0.008 


0.02 


Look Angle 


0.001 


0.02 


0.003 


0.02 


Sector 


0.01 


0.02 


0.005 


0.02 



Note. — Given are fractional errors relative to the total sky brightness S = S p + S s . 



P10 S% lpp fluxes, i.e., S Rlpp = 2S£ Ipp for the P10 mea- 
surements. We group each exposure with the neigh- 
bors within 5°, and reject it if it deviates more than 
three times the standard deviation from the median of 
the group either in Sb ipp , Sr ipp , or Sb ipp / Sr ipp ■ This 
process reveals that the P10 measurements on the day 
216/1973 are very noisy, with nearly 50 % rejected, in 
contrast to at most ~10% rejected for other days. There- 
fore all the data taken on this day are discarded. Table 
Q] lists fractions of the usable exposures as well as the 
total number of exposures taken on each day. We find 
that approximately 70% of the IPP data are usable for 
the present purpose. 

2.3. Mapping Process 

We integrate all the usable data into a single sky map 
at each band. Since IPP FOVs follow complex trajec- 
tories on the celestial sphere and most of the sky were 
observed several times or more, the final maps have much 
higher sensitivity and spatial resolution than the individ- 
ual measurements. The map generation is limited to the 
high Galactic latitude sky (|6| > 35°) where foreground 
emissions from Galactic sources are relatively faint. 

First, a set of pixels is placed on the celestial 
sphere based on the Galactic coordinate system, us- 
ing Lambert's zenithal e qual-area projection (e.g., 
ICalabretta fc Greisenl I2002D . I t is the same projection 
as used bv lSchlegel et al.l ( )1998l ) for creating the 100 ^m 
emission map described below. Pixel size of 0°.32 x 0°.32 
is adopted, which projects all-sky data onto two 512 x 
512 pixel maps centered on the north and south Galactic 
poles. The brightness Sij at each pixel (i, j) is calculated 



by 

Si* = <S&>, (1) 

where Sfj represents an ensemble of the n = 1, 2, iVth 
IPP measurements whose FOVs contain the pixel (i, j) 
and (S™j) is mean brightness of the ensemble. When cal- 
culating (S™j), we apply the following rejection criteria 

based on diffuse emission brightness S dlSuse = S — 5 star , 
where S stliI is starli ght brightness in each FOV calcu- 
lated in Section [330. We reject the data deviating more 
than three times the standard deviation from the mean of 
the ensemble either in S£ iffusc , S'^ iffusc , or Si 1 ^ / St iSusc . 

dipp ' Jiipp ' dipp ' nipp 

The unreliable measurements with S < S stm ' — 20 bgu 
are also rejected (see Section |3~T|) . Associated error at 
each pixel is estimated from the scatter of the ensem- 
ble Sfy This mapping algorithm is similar to first it- 
eration of that adopted by I Gordon et al.l (|1998f ). and 
to first iteration of th e "maximum correlation method" 
(jAumann et al.lll990l ) used for improving angular reso- 
lution of the Infrared Astronomy Satellite (IRAS) maps. 
Considering relatively low signal-to-noise ratios at low 
brightness regions where our study focuses on, we do not 
perform the higher-order iteration of the map generation 
algorithm in the present analysis. 

We show the created -Brpp- and i?ipp-band maps of the 
high Galactic latitude sky (|6| > 35°) in Figures [3] and [U 
The above process projects ^40,000 IPP exposures onto 

4 It is better to use 5 dlf T use for the rejection process than to 
use S itself because the latter can deviate significantly from each 
other in the ensemble S™. due to different bright stars falling on 
the FOVs. 
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-200,000 pixels (-20,000 deg 2 ) at the high Galactic lat- 
itudes, which gives the rough estimate of the resultant 
angular resolution — 0°.7. On average, 16 measurements 
contribute to each pixel except in the fields with no ob- 
servations. Finally we apply 4x4 pixel binning to the 
created maps. 

At this stage we can further check accuracy of the IPP 
measurements by investigating reproducibility of the cre- 
ated maps. We divide the data into two subsets based 
on each of the following five criteria, which are selected 
to maximize systematic difference between the subsets; 
spacecrafts (P10 or Pll); the heliocentric distances R 
(< 4.5 AU or > 4.5 AU); elongations (< 115° or > 115°); 
look angles (\L - 90°| < 30° or > 30°); and sectors (01 
- 32 or 33 - 64). The demarcations of the heliocentric 
distances (4.5 AU), elongations (115°), and look angles 
(30°) are chosen so that the divided two subsets have 
similar numbers of measurements. Then a pair of sky 
maps are created adopting each of the five divisions, and 
distributions of the brightness deviation 

cdivl cdiv2 cdivl c cdiv2 c 

J iJf O iD iD lJ 

= 5 divl + 5 div2 = S = S 

are calculated in the overlapping regions. Here, g dlvl 
and S dlv2 are brightness of same pixels of the two maps 
and S = (S" divl + S div2 )/2. The comparison is limited to 
the pixels for which more than three measurements con- 
tribute and within the brightness range S — 30 - 90 bgu. 
The results are shown in Figure [5j All the distributions 
are close to the Gaussian distribution; the resultant esti- 
mates of systematic and random errors from the best-fit 
mean and standard deviation values are summarized in 
Table [2] The table shows that any subset of the data 
reproduces the final sky map within 1 % systematic un- 
certainty in both of -Bipp and i?ipp. Thus, we conclude 
that systematic errors inherent in the IPP measurements 
are within 1 % of the total sky brightness. 

3. GALACTIC LIGHT 
3.1. Starlight 

Outside the detectable ZL clouds, the dominant 
brightness component incident on the IPPs is Galactic 
starlight. Contribution of bright stars has already been 
subtracted when the Background Sky Tape was created, 
referring to stars in the Yale Bright Star Catalog and 
V < 8 stars in the USNO Photoelectric Catalog (see 
Section f2 . 1 [) . Therefore, we have to subtract contribu- 
tion of fainter stars from the IPP measurements. For 
this purpose, we use two all-sky star catalogs along with 
a star-count model as described below. Transformations 
of stellar colors between the different filter systems are 
also developed. 

Integrated brightness of relatively bright (V < 11 
mag) stars are c alculated using the Tycho-2 Catalog 
(|H0g et al.ll2000aD . The catalog is based on observations 
by the star mapper of the ESA Hipparcos satellite, and 
contains positions, proper motions, and photometry data 
for —2,500,000 brightest stars in the sky. Source magni- 
tudes were measured at the two passbands Bt and Vt, 
which are close to the Johnson B and V bands. Detection 
completeness is —90 % at V = 11.5 mag, with magnitude 
errors less than 0.1 mag for brighter stars. Systematic 
uncertainty of the photometry is estimated to be smaller 



than 0.02 mag (|H0g et al.ll2000bD . 

Contribution of fainter (V > 11 mag) stars is derived 
from the HST Guide Star Catalog II (GSC-II) version 
2.3 (|Lasker et all 120081 ). It was constructed from the 
scanned images of 9541 Palomar and UK Schmidt pho- 
tographic sky survey plates digitized at Space Telescope 
Science Institute (STScI). About 950,000,000 sources to 
the limiting magnitude of R — 20 mag are contained in 
the catalog with photometry information mainly at the 
three natural plate passbands Bj, Rp, and In- The lim- 
iting magnitudes of the catalog correspond to -Bipp = 
22.5 mag and i?ipp =20.5 mag (see below for transfor- 
mations of stellar magnitudes). Magnitude errors are less 
than 0.2 mag at Rf < 19.5 mag, and systematic uncer- 
tainty of the photometry is estimated to be within a few 
hundredths of magnitude. 

We estimate contribution of stars even fainter than 
the GSC-II detection limits using a star-count model 
provid ed by a stellar popu lation synthesis code TPJLE- 
GAL (IGirardi et al.ll2005h . TPJLEGAL generates and 
spatially distributes stars according to the user-defined 
parameters, such as star formation rate, age-metallicity 
relation, initial mass function, and distribution geome- 
try, and then predicts stellar photometry at any field 
of the Galaxy (also at star clusters and nearby galax- 
ies). We adopt the default set of parameters for deriving 
Galactic star counts at the Bessel B and R bands down 
to 32.0 mag. Then the ratios S iaint / S hright are calcu- 
lated at 30 representative Galactic coordinates, where 
gbnght anc j gfamt re p resen t integrated brightness of stars 
in the brightness ranges -Bipp (i?ipp) = 15.0 - 22.5 (15.0 
- 20.5) mag and 22.5 - 32.0 (20.5 - 32.0) mag, respec- 
tively. Note that the demarcations between bright and 
faint stars correspond to the GSC-II limiting magnitudes. 
We find that these faint star fractions are small, rang- 
ing from Shinty ^bright = q.01 to 0.03 at Bipp and from 
^faint^bright _ q Q3 to 0.12 at R 1PP . The ratios vary 
fairly smoothly between the 30 positions, hence we cal- 
culate those at any positions by interpolating among the 
three neighbors chosen from the 30 points. The effects 
of the uncertainty in the above process for estimating 
contribution of the faintest stars on our final results are 
discussed in Section 5. 

Transformations of stellar magnitudes from the Hip- 
parcos (Bt and Vt), the natural plate band (Bj, Rp, 
and In), and the Bessel (B and R) systems to the 
IPP (-Bipp and i?ipp) system are developed using the 
Bruzual-Persson-Gunn-Stryker (BPGS) stellar spectral 
atlas provided by STScI. The atlas contains the optical- 
to-near-IR spectra of 175 stars with a wide range of spec- 
tral type and luminosity class, and is often used to simu- 
late color-color dia grams of various Galactic stars (e.g., 
iHewett et al.l l2~006). Figure [6] shows the relations of stel- 
lar colors in the different passband systems. All the re- 
lations can be well described by one- or two-component 
regression lines, which are used for converting stellar col- 
ors in this work. For convenience, we summarize the 
coefficients of the regression lines in Table [3] 

Using the above star catalogs and star-count model 
with the stellar color conversions, we compute integrated 
brightness of starlight S star falling on FOV of each IPP 
measurement. In order to avoid duplicated counting, all 
the stars found in the Yale Bright Star Catalog and the 
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Fig. 3. — Pioneer 10/11 IPP Bipp-band maps of the north (left) and south (right) high Galactic latitude sky (|ft| > 35°). Color is linearly 
scaled between 30 (black) and 75 (white) bgu. The black lines show the Galactic longitudes of / = 0°, 90°, 180°, and 270°, while the black 
circles show the Galactic latitudes of b = ±60° . The Galactic poles are at the centers of the two images. 



North (b > 35) 
270 



South (b < -35) 
90 



180 



180 



| 



1 = 90 



1 = 270 



Fig. 4. — Same as Figure[3] but for the IPP Hipp band. Color is linearly scaled between 35 (black) and 80 (white) bgu. 



V < 8 stars in the USNO Photoelectric Catalog are 
excluded from the Tycho-2 and GSC-II catalogs. The 
Tycho-2 stars in the GSC-II are also removed. Fractions 
of the fainter stars than the GSC-II detection limits de- 
rived from the TRILEGAL code, S iaint / S hTisH , are used 
in combination with the actual values of 1 5 ,brlght to es- 
timate contributions of the faintest stars down to 32.0 
mag. Then the starlight brightness is subtracted from 
each IPP FOV, resulting in diffuse emission brightness 
^diffuse _ g _ gsta.i_ Random uncertainties in the above 
process of starlight estimation and subtraction, including 
those originating from the relatively large IPP pointing 
uncertainty, are directly reflected in the final data scatter 
as we see in Section 4. 
At this stage, we find that a small fraction (^5 %) 



of the IPP brightness are significantly lower than the 
starlight brightness of the fields, i.e., 5* < 5' star — 20 bgu, 
which are roughly more than three times the standard 
deviation of the data scatter. Since there are few data 
showing the opposite deviation S > S staY + 20 bgu, it is 
most likely caused by some failures in the IPP measure- 
ments rather than the natural data scatter. These data 
are rejected in the map generation process. 

The diffuse emission brightness of IPP FOVs are in- 
tegrated into sky maps using the same algorithm as de- 
scribed in Section [2~3l i.e., 



5" 



diffuse 



diffuse;n\ 
3 ' 



(2) 



where S', 



diffuse;', 
i,3 



is an ensemble of the n — l,2,...,7Vth 



IPP diffuse emission brightness whose FOVs contain the 
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TABLE 3 

Conversion of Stellar Colors 



Color 1 (To) 


Color (From) 


a 


b 








Filter System (From) 




-R F 


Bj 


-R F 


0.00 


1.18 








GSC-11 north 


Bj- 


Ripp 


Bj 


- R F 


0.00 


0.96 








GSC-II north 


B\pp 


-In 


Bj 


-In 


0.03 


1.11 








GSC-II north 


Bj- 


Ripp 


Bj 


-In 


0.01 


0.72 


at 


color 


< 2.66 


GSC-II north 










1.23 


0.26 


at 


color 


> 2.66 


GSC-II north 


B IPP 


-In 


Hp 


-In 


0.07 


4.56 


at 


color 


< 0.64 


GSC-II north 










2.09 


1.41 


at 


color 


> 0.64 


GSC-II north 


■Ripp 


-In 


Rp 


-In 


0.03 


1.04 








GSC-II north 


Bipp 


- R F 


Bj 


- R F 


0.00 


1.23 








GSC-11 south 


Bj- 


Ripp 


Bj 


— Rp 


-0.01 


1.00 








GSC-II south 


Bjpp 


-In 


Bj 


-In 


0.04 


1.13 








GSC-II south 


Bj- 


-Ripp 


Bj 


-In 


0.02 


0.71 


at 


color 


< 2.54 


GSC-II south 










1.11 


0.28 


at 


color 


> 2.54 


GSC-II south 


Bppp 


-In 


Rp 


-In 


0.05 


4.23 


at 


color 


< 0.68 


GSC-II south 










2.00 


1.38 


at 


color 


> 0.68 


GSC-II south 


Ripp 


-In 


Hp 


-In 


0.01 


0.98 








GSC-II south 


Bipp 


- V T 


B T 


- V T 


0.03 


0.85 








Tycho-2 


Bjp — 


-Ripp 


13 p 


- V T 


-0.02 


1.64 








Tycho-2 


B\pp 


-Br 


B 


-R 


0.02 


1.00 








Besscl 


Bg - 


Ripp 


B 


-R 


-0.03 


1.04 








Bessel 



Note. — Conversion from the color to the color 1: [color 1] = a + b x [color 0]. 
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Fig. 5. — Distributions of the brightness deviation D = (S dlvl — 
S div2 )/(5 divl +S div2 ) at the B IPP (panels (a)-(e)) and R IPP - (pan- 
els (f)-(j)) bands (histograms). Division of the data are based on 
spacecrafts (panels (a) and (f)), the heliocentric distances (panels 
(b) and (g)), elongations (panels (c) and (h)), look angles (panels 
(d) and (i)), and sectors (panels (e) and (j)). See the text for the 
full description. The dotted lines are drawn at D = 0. The super- 
posed smooth curves are the Gaussian distributions fitted to the 
data. 

pixel (i, j). The same rejection criteria as used in Equa- 
tion ((T|) are applied, and associated error at each pixel is 



estimated from the scatter of the ensemble S', 



diffuse;', 
i,3 



Fi- 



nally, 4x4 pixel binning is applied to the created diffuse 
emission maps. 



3.2. Diffuse Galactic Light 
3.2.1. Overview 

The DGL is attributed to dust and gas particles in 
the interstellar medium (ISM) . Dominant contribution to 
the optical DGL comes from scattering of the interstel- 
lar radiation field (ISRF) by the dust. This scattering 
process is the p rime cause of the well-k nown interstel- 
lar ex tinction. iBernstein et al.l ()2002aD and iBernsteinl 
(2007) derive brightness of this scattered DGL compo- 
nent in their field (Galactic coordinate I = 206°. 6 and 
b = —59°. 8) by combining a simple scattering model and 
thermal emission brightness at 100 /im wavelength. The 
estimated value is ~ 1 bgu, which is a few tenths of the 
COB brightness claimed by them. 

Another process generating the DGL is line and con- 
tinuum emissions from warm (T ~ 10 4 K) ionized gas in 
the ISM. The strongest line emission contributing to the 
optical DGL is Ha, which falls in our i?ipp band. The 
fainter emission lines [N II] A6583 and [S II] A6716 could 
also contribute to th e DGL, each wit h roughly ~0.3 times 
the Ha brightness (|Revnoldslll985D . Hydrogen Balmcr 
lines other than Ha (H/3, H7, ...) could contribute to 
the DGL in the -Bipp band. Ado pting the e stima tes of 
diffuse Ha intensity provided bv iRevnoldsl (|1992f ). the 
total DGL brightness due to the diffuse line emissions 
at |6| = 60° are less than 0.2 and 0.3 bgu at -Bipp and 
i?ipp, respectively (the case B condition of the recombi- 
nation theory is assumed). On the other hand, contin- 
uum emissions from warm ionized gas such as the two- 
photon, fre e-free, and boun d -free e missions, are almost 
negligible. IBernstein et al.l (|2002af l estimate the DGL 
brightness due to the continuum emissions as less than 
0.01 bgu at the optical wavelengths longer than 0.37 /xm. 

In total, the DGL brightness at -Bipp and i?ipp are 
at most a few bgu at \b\ ~ 60° and less at the higher 
Galactic latitudes, most of which comes from scattering 
of the ISRF by the interstellar dust. 

We estimate the DGL contribution to IPP diffuse emis- 
sion brightness by separating out the component which 
correlates with the diffuse Galactic far-IR emission. It is 
well known that the far-IR emission traces thermal ra- 
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Fig. 6. — Conversions of stellar colors between the Pioneer IPP and the star catalog systems: the Palomar natural plate passbands for 
the GSC-II northern stars (panels (a)— (f)), the UK Schmidt natural plate passbands for the GSC-II southern stars (panels (g)— (1)), the 
Hipparcos passbands for the Tycho-2 stars (panels (m) and (n)), and the Bessel passbands (panels (o) and (p)). The open circles represent 
the 175 BPGS stars while the dashed lines show the regression lines. 



diation from the interstellar dust heated by the ISRF. 
Hence it is brighter in the fields with larger column den- 
sity of the dust (and gas) and stronger ISRF, where the 
DGL is also brighter. It naturally results in correlation 
between the DGL and far-IR brightness especially in op- 
tically thin dust regions, where the both brightnesses are 
proportional to the dust column density. Actually sev- 
eral authors find linear correlations between the DGL 
and t he diffuse Galactic 100 ulvh brightness in different 
fields (iLaureiis et al.lll987tlGuhathakurta fc Tvsonlll989t 
iPalev et all 119911 ). while the reported DGL-to-100 /mi 
brightness ratios vary from field to field. As well as non- 
uniform dust properties, large measurement errors inher- 
ent in the treatments of the foreground emissions are ob- 
viously one of the major contributors to this field-to-field 
variation. With the DGL-to-100 /im brightness ratio a<j, 



the IPP diffuse emission brightness can be written as: 



^diffuse _ gDGL _|_ gCOB 
— fldWloO/jm 



diffuse qCIB \ , qCOB 



(3) 



where S DGh and S COB are the DGL and COB brightness 
at the IPP bands, and S($™£ is diffuse 100 /im bright- 
ness observed outside the ZL clouds. The cosmic infrared 
background (CIB) component Sfrx^m is subtracted from 

•^lOOurn m the above equation, leaving only Galactic com- 
ponent. 

3.2.2. Simulation 

In order to make a proper interpretation of the ob- 
served relation between g^" 80 and S'foo"^ in the pres- 
ence of the data scatters, here we perform a simple Monte 
Carlo simulation. We adopt the numbers and distribu- 
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tions as similar as possible to what is observed (found in 
Section 4) , while we do not intend to deduce any accurate 
numbers from this simulation. The COB and CIB bright- 
ness are assumed to be S COB = 1.0 bgu and S^qq^ = 0.8 

MJy sr -1 . We generate ^8000 mock data with 100 /im 
brightness following the simple probability distribution ex 

av^(— cdiffuse ^ Q+ odiffusc ^ gCIB „rV,;„V, voo,,l+c in 



which results in the 



exp( -SVqq^) at "SVocvm ^ '-'loo/iim 
similar number distributions of S^qq^ to the observed 
at the high Galactic latitudes. 5 fdlffusc is calculated by 
Equation (j3|) assuming the Gaussian probability distri- 
bution of ad with the mean a d and the standard deviation 



, ad , where a d 



3.3 bgu (MJy sr 



and a ad =1.0 bgu 



(MJy sr -1 ) -1 are adopted. The former value is similar 
to what is found in the actual measurements while the 
latter is coupled with measurement errors and is hard to 
derive from the present data. The resultant input rela- 
tion between 5" dlffuse and S'fgQ^j is shown in the top panel 
of Figure [Jj Then the Gaussian random errors are added 
to S diSusc and Sffloj^ with the standard deviations of 5.0 



bgu and 0.2 MJy sr 
similar S diSuse " "" 



odiffuse 
°100^m 



respectively, which reproduces a 
plot to the observed. We calcu- 
late mean values of 5' dlffusc and their errors in Sfoo^m bins 
stepped by 0.1 MJy sr -1 , which are shown in the bottom 
panel of Figure [Jj (red circles) . We repeat the simulation 
100 times varying the random components, and obtain 
mean <S' dlfTuse -S' d Qo^ relation of the whole realizations as 
plotted in the same panel (green diamonds). 

Figure [JJ shows that in the presence of the realistic 
data scatters (red circles), we will find the two dis- 
tinct features in the observed S diffusc ~S d k^ plot. At 



diffuse 
100/rni 



> 



cCIB 
a 10Q>m> 



_ cdittusc 
°lu0/xm 

the original correlation of the DGL 
component is simply reproduced with some random de- 
viations. This cor relation breaks at Sfffij^ ~ S^ m , 
and a relatively flat distribution of S dlSuse is seen at 
the smaller Sfjj$£. The S ditlusc value in this region is 
dominated by the COB brightness while a minor con- 
tribution of the DGL is present. Hence, such a inflec- 
tion point of the 5 ldlffuse -5' d oo)im relation gives rough es- 
timates of the COB and CIB brightness. Alternatively, if 
the CIB brightness S , p I |^ lm is known in advance, we can 
obtain more accurate estimates of the COB brightness at 
the intersect of the regression line of the correlation at 



cdiffuse ^ cCIB 
o 100>m ^ J 100Atm 



and the line of S$$£ 



^loc^m- Note 



that it corresponds to substituting Sf^^ 
Equation ©. 

4. RESULTS 

We show the measured IPP diffuse emission bright- 
ness ,S dlffuse versus the diffuse 100 jum brightness S d oo)}m 
outside the ZL clouds in Figure [5] The present analysis 
focuses on the lowest brightness region with Sfoopm < 3.0 
MJy sr -1 and with the IPP map coverage at the Galactic 
latitudes |6| > 35°, which corresponds to about a quarter 
of the whole sky. Mean values of 5' dlffusc and their errors 
are calculated in S^o"^ bins stepped by 0.1 MJy sr -1 
(red circles), which are used for the regression analysis in 
what follows. The errors are derived from the data scat- 
ters, which arise from the following three main causes: 
the IPP measurement errors, uncertainties in the esti- 
mation and the subtraction of Galactic starlight, and in- 
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Fig. 7. — Input (panel (a)) and output (panel (b)) S diffusc and 
•SioOfim of one realization of the Monte Carlo simulation (dots). 
The dotted lines show the assumed COB and CIB brightness, while 
the dashed lines show the mean DGL-to-100/im relation adopted 
(Equation J3} with a d = a d ). The red circles and error bars rep- 
resent mean values of the output S dlSuse and their errors in the 
■^loojlm Dms of the realization, while the green diamonds represent 
mean values of all the 100 realizations. The sizes of the red circles 

^diffuse 
'lOO/jm 



trinsic scatters of DGL-to-100 /im brightn ess ratios. The 
diffuse 100 /im brightness is taken from Schlcgcl ct al. 
(fl998h . They reprocess and combine the IRAS and the 
Cosmic Background Explorer (COBE) Diffuse Infrared 
Background Experiment (DIRBE) measurements which 
results in an all-sky 100 /im brightness map with high 
angular resolution provided by the IRAS and accurate 
absolute flux calibration provided by the DIRBE. The 
ZL and confirmed point sources and galaxies have been 
removed, while no significant CIB is detected in their 
analysis. We apply 32 x 32 pixel binning to their 100 
/im map so that it has the same resolution as our IPP 
maps. Uncertai nty of the 100 /im b rightness is not ex- 
plicitly given in iSchlegel et a l. (1998), while accuracy of 
the interstellar dust reddening (Eb-v) derived from the 
100 /im map is estimated to be 16%. Since Eb-v has 
a larger uncertainty than the original 100 /im brightness 
due to the various assumptions made, e.g., those made 
in estimating dust temperatures, accuracy of the 100 /im 
brightness would be considerably better. In addition, the 
32 x 32 pixel binning applied in the present work should 
further reduce random errors. 
While the better 100 /im brightness map is pre- 
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sented bv lMiville-De schcnes & Lagache (2005), their Im- 
proved Reprocessing of the IRAS Survey (IRIS) map 
includes point sources and thus cannot be compared 
directly to our diffuse emission maps. However, a 
very good ag r eemen t is found between the IRIS and 
ISchlegel et al.l (119981) data sets at scales larger than 
~ 0° 5 (|Miville-Deschenes fc Laeachd I2005D where the 
present analysis focuses on (note the IPP map resolution 
of ~ 0°.7). Hence, we do not expect a significant change 
of the final results when the IRIS map is alternatively 
used after point-source removal. 

In Figure[51 we clearly detect the linear correlations be- 
tween S diSnsc and Sf^Z at the large Sf$$Z, and more 
importantly, the expected breaks at S d ^ s £ < 0.8 MJy 
sr -1 at the both IPP bands. These inflection points are 
in very good agreement w ith the CIB brightness reported 
by ILagache et all (l2000h . Sf ™m = 0.78 ± 0.21 MJy 
sr , who made use of the same COBE/DIRBF, mea- 
surements as used for absolute calibration of the present 
100 fim map. However, we note that S^^^ here can in- 
clude any brightness component which does n ot correlate 
with t he interstellar dust emission. Actually IDole et al.l 
(12006T) point out tha t the CIB brightness measured by 
ILaeache et al.l (|2000f) is likely to include the residual ZL 
01 S wo£n ~ 03 MJ y sr_1 - The present results indicate 
that a similar a mount of the residual ZL is present in the 
100 iim map of ISchlegel et all (fl99ll . 

Our estimates of the COB brightness are derived as 
the S diSusc values on the observed 5 diffuse -5^f^ linear 
correlations made by the DGL, at the point where 100 



/zm brightness equals to the CIB (i.e., S 



diffuse 
lOO^m 



CIB 
100/jm 



in Equation ([3])). First, we obtain the regression lines 
of the 5 diffusc -S'^ s n ° relations at Sf$,™£ = 1.0 - 3.0 
MJy sr -1 using the least-% 2 method. The derived mean 
DGL-to-100 fim brightness ratios are a° d = 3.2 ± 0.1 and 
3.4 ± 0.1 bgu (MJy sr -1 ) -1 at £?ipp and i?ipp, respec- 
tively. Then the intersects of the regression lines and 
Sfoo^m = S iao»m are derived, whose S diSusc values give 
our COB es timates. Here, we adopt the CIB brightness 
reported bv lLagache eTaLl (pOOOl ) (S^ m = 0.78 ± 0.21 
MJy sr -1 ) and include the associated uncertainty in the 
s cob errors . The reS ults are S COB = 1.8 ± 0.7 and 1.2 
±0.7 bgu at Bipp and j ?ipp, respectively. 

Recently, IPenin et al.1 (j 2 1 lh present a CIB brightness 
01 S ?wvm = °-20 ± 0.09 MJy sr" 1 from the analysis of 
the Spitzer Space Telescope (Spitzer) Multiband Imag- 
ing Photometer for Spitzer (MIPS) data. Based on the 
MIPS calibration at 160 ^im, the 100 /jm CIB is de- 
rived from the CIB anisotropy color between 100 and 160 
^m computed from the MIPS and the IRIS data. The 
reason for t he dis crepancy betwee n their result and the 
IDole et all (|2006ft revision of the iLagache et all (|2000l ) 
result, 0.48 ±0.21 MJy sr -1 , remains unclear. However, 
the latter is the better choice for the present study, since 
it is based on the same DIRBE calibra tion as used for 
the 100 (im map of ISchlegel et al.l £[998). For reference, 
the estimated COB brig htness would be red uced by ^0.9 
bgu when the result of IPenin et al.l (|2011f ) is used with 
the plausible residual ZL level of S{ e f^ m - 0.3 MJy sr" 1 
in the 100 fim map (see above). 




-DIRBE 




IRAS-DIRBE 

Fig. 8. — Observed IPP diffuse emission brightness S d 



at 



Bipp (panel (a)) and -Ripp (panel (b)) vs. the diffuse 100 (im 
brightness Sioo^m (^°* s ) ■ '^^ le re< ^ circles and error bars represent 



mean values of S 



and their errors in the S? 



bins. The 



J 100Mm 

sizes of the circles are proportional to the numbers of the data 
points in the bins. The solid green lines show the regression lines 
at Sf$*%i > 1.0 MJy ST' 1 , while the dashed green lines show the 
CIB brightness 5f : ( f Mm reported by ILagache et al.l l l2000Tj . The 
shaded areas show 1<t confidence intervals of our final estimates of 
the COB brightness. 

The above 5" COB errors do not include associated sys- 
tematic uncertainties. Median IPP brightness of the 
fields with Sf^j^ = 0.78 ± 0.21 MJy sr" 1 before starlight 
subtraction is S ~ 40 bgu at the both bands, which leads 
to additional error of 0.4 bgu in the presence of 1% sys- 
tematic uncertainty in the IPP measurements as esti- 
mated above. The Tycho-2 and GSC-II stars contribute 
about equally (i.e., ~ 20 bgu each) to the above starlight 
while fainter stars make negligible contribution (see Sec- 
tion 5). Taking 0.02 mag systematic uncertainties of the 
stellar magnitudes into account, 0.5 bgu error is further 
added. After combining all these errors, the final esti- 
mates of the COB brightness are S COB = 1.8 ± 0.9 and 
1.2 ± 0.9 bgu at -Bipp and i?ipp, respectively, or 7.9 ± 4.0 
and 7.7 ± 5.8 nW m~ 2 sr -1 . The results are summarized 
along with the derived a° d values in Table |U 

5. DISCUSSION 

5.1. Correlation between the DGL and the 100 [im 
Emission 

By converting bgu unit to MJy sr -1 , we obtain di- 
mensionless representations of the mean DGL-to-100 /im 
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TABLE 4 

COB Brightness and Mean DGL-to-100 /im Brightness Ratios aj 



Wavelength COB Brightness DGL-to-100 fim Ratio a d 

Band (/-im) (bgu) (nW m — 2 sr _1 ) (bgu [MJy sr — 1 ] _1 ) (dimensionless) 

Bipp 044 1.8 ± 0.9 7.9 ± 4.0 3.2 ± 0.1 (2.1 ± 0.1) xl0~ d 

giPP 0-64 1.2 ± 0.9 7.7 ± 5.8 3.4 ± 0.1 (4.6 ± 0.1) xl0~ 3 



brightness ratios a° d = (2.1 ±0.1) x 1(T 3 and (4.6 ±0.1) x 
10~ 3 at i?ipp and ifjpp, respectively, which are useful for 
a comparison with other studies. These values are also 
found in Ta ble IH The previ o us mea surements of a d are 
compiled by lBernstein et ahl (|2002al ) in their Table 7. 

At the blue b and around 0.45 /mi, 
IGuhathakurta fc Tyson! ([1989) find the values of 
a d = 1.1 x 10" 3 and 2.6 x 10" 3 in their low 100 /mi 
brightness regions "ir 2" and "ir 3". They can be 
reproduced by a simple model of the interstellar dust 
scattering the ISRF, which predicts a d — 3.3 x 10~ 3 . 
A similar value of = 1.1 x 10~ 3 is observed by 
iPalev et al.1 (|1991[ ). IBernstein et~aLl (|2002al ) provide 
another dust scattering model and show that a d — 
1.6 x 10~ 3 at the high Galactic latitudes. From the ob- 
servations of optic ally thin high Galactic latitude clouds, 
iWitt et alJ (|2008l) find the ratio of a° d = 2.1 x 10~ 3 . 
Overall, our result at £?ipp agrees with these previous 
measurements and model predictions. 

At the red band around 0.65 /im, the measured a d val- 
ues are 2. 2 x 10~ 3 and 4.4 x 10~ 3 in the i r 2 and ir 3 
region s of IGuhathakurta fc Tvsonl (|1989f ). IPalev et all 
CEMD) report a° d = 8.5 x 10" 3 at 0.55 /im and 11 x 10" 3 
at 0.70 /tm. These measurements point the a d (R)/a d (B) 
ratios of at least 2.0, while they are known to be sig- 
nificantly larger than the predi ctions of dust scattering 
model s. In fact the model s of IGuhathakurta fc Tvsonl 
(1989) and IBernstein et all (|2002a|) consistently predict 
a d (R)/a d (B) ~ 1.4. The excess component in the 
red band is called extended red emission (ERE) and 
is found in many Galactic objects including refle ction 
nebulae (IWitt et all 119841: IWitt fc Borosonl Il990l). in- 
frared cirrus (IGuhathakurta fc Tyson 1198^7 . and plan- 
etary nebulae jFurton fe Wittl 119901 Il992h . Its origin 
is most likely the photoluminescence process of mate- 
rials such as hydrogenated amorphous carbon and poly- 
cyclic aro matic hydrocarbon excited by UV/optical pho- 
tons (see IWitt fc Viihl 12004 for a review). Our result, 
a d (R)/a d (B) — 2.2, is consistent with the previous ob- 
servations, suggesting that ERE is also present in the 
diffuse ISM with t he lowest far-IR brig htness. It con- 
firms the finding of iGordon et al.1 (|1998f) who reach the 
same conclusion from the analysis of the IPP data. 

In summary, our results are in overall agreement with 
the previous observations toward the denser dust regions. 
A further study on this issue is beyond the scope of this 
work and will be presented in a future paper. 

5.2. Resolved Fraction of Cosmic Background 

We compile the current measurements of the cosmic 
background and the integrated brightness of galaxies 
at ultraviolet (UV), optical, and near-IR wavelengths 
in Figure [9] At UV wavelengths, the upper limits 
of the cosmic background are obtained from the anal- 



ysis o f the Voyager Ultra violet Spectrometer (UVS) 
data (jEdelstein et al.l I2000D and from the H5T/Space 
Telescope Imaging Spectrograph (STIS) observations 
(|Brown et al.l [2000). They are a few times the inte- 
grated brightness of g alaxies measured by the HST/STIS 
(jGardner et al.l 120001 ) . still leaving a large gap to be 
bridged. The situation i n near-IR wavelength ra nge is 
much more controversial. iMatsumoto et ahl ((2005! ) claim 
detection of the strong near-IR CIB based on the In- 
frared Telescope in Space (IRTS) data. Their CIB val- 
ues are marginally consistent with the results from the 
C O-Rfi yPIRBE measurements reported by several au- 
thors (IGoriian et al.ll2000t IWrightJl200U ICambresv et al.l 
120011: 1 Wright! I2004D at 1.25, 2.2, and 3.3 /tm. The in- 
tegrated brightness of galaxies at these wavelengths, as 
well as at the optical UBVI bands, a re derived from the 
Hubbl e Deep Field (HDF) data set bv lMadau fc Pozzettil 
( 2000) . Those at the four bands of the Spit zer Infrared 
Array Camera (IRAC) are presented by Fazi o et al.1 
(2004). Usi ng the HDF and t he Subaru Deep Field 
(SDF) data iTotani et al.1 (120011) obta in the consistent 
results with iMadau fc Pozzettil (|2000D . They also find 
that 80%-90% of the total light from normal galaxies has 
already been resolved in the SDF J and K bands, based 
on a galaxy evolution model taking into account vari- 
ous selection effects of observations. Therefor e, the large 
CIB excess found by IMatsumoto et all (|2005|) should be 
attributed to either some exotic radiation sources such 
as Population III stars or the residual ZL in the IRTS 
data. Another constraint comes from the High Energy 
Stereoscopic System (HES S) 7-ray observa t ions of the 
blazars at z = 0.17-0.19 bv lAharonian et ahl (|2006l) . who 
find that the unexpectedly hard spectra of these objects 
cannot be reproduced without assuming very low upper 
limits of the CIB. Their upper limits are close to the in- 
tegrated brightness of galaxies and are in sharp contrast 
with the previous CIB est imates. 

At optical wavelengths, IBernstein et al.l (|2002aD claim 
the "first detections" of the COB at 0.30, 0.55, and 0.80 
/im (see Section 1). They used the Wide Field Plane- 
tary Camera 2 (WFPC2) on board the HST for abso- 
lute surface brightness measurements of the night sky, 
and a ground-based telescope for contemporaneous mea- 
surements of the ZL b rightness . How ever, a number of 
problems are raised by iMattilal ([2003) for their analysis, 
result ing in a series of corrections of the p resented re- 
sults (IBernstein et al.1 120051: IBernstein! 12007! ) . Their lat- 
est values are plotted in Figure [9j which are nearly one 
order of magnitude larger than the integrated brightness 
of galaxies and seem to be hard to explain in the current 
framework of galaxy evolution. 

Our new results of the COB brightness are mildly 
larger than and consistent within la uncertainty with 
the integrated brightness of galaxies. They are on 
the smooth extension of the upper limits found by 
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Fig. 9. — Current measurements of the cosmic background (filled symbols) and the integrated brightness of galaxies (open symbols) 
at UV, optical, and near-IR wave lengths. The cosmic background measurements include the UV upper limits (blue arrows) at 0.10 (im 
obtained from the Voyager/WS (Edclstcin ct al. 2000) and at 0.16 fxra from the HST/STIS (Brown ct al. 2000), the claimed detections 
at optical wav e length s using the HST/WFPC2 (Bernstein 2007j, green s quares) and at n e ar-IR wavelengths usin g the C O-RE/DIRBE 
|Goriian et"aTI (2000), green diamonds; IWrighH (120011) . purple diamonds; Cambrcsv et al.l i 20011) . blue diamonds; Wright (2004 j), grey 
diam onds; the wavelengths of these measurements are slightly shifted relative to each other for clarity] and the IRTS (Matsurnot o et al.1 
2005, black circles). The red star s are the Pioneer/TPP results of this work, while the red solid line with arrows between 0.8 and 4 /im 
represe nts the HESS upper limits l lAh aronian ct al. 2006). The integrated brig htness of galaxies come fr om the HST/STIS measurements 
at UV (Gardner ct al. 2000, squares), the HDF compilation from UV to near-lR (Madau & Pozzctti 2000, triangles), and the Spitzer/IRAC 
measurements at near-IR wavelengths (Fazio ct al. 2004, diamonds). 



lAharonian et al.l ((2006). With our best estimates, ap- 
proximately 60% and 90% of the COB have already been 
resolved into discrete galaxies in the HDF at 0.44 /im 
(.Btpp ) and 0.64 (im (Rrpy ), respectively. On the other 
hand, iTotani et all (|200lf ) demonstrate that 60%-90% 
and 80%-100% of the total light from galaxies have been 
resolved at 0.45 and 0.61 /im. The above facts indicate 
that bulk of the COB are comprised of normal galaxies, 
and there are little room for contributions of other pop- 
ulations at these wavelengths. Note that AGNs are one 
of such populations, while their contributions to the cos- 
mic background are estim ated to be less than 10%-20% 
(jMadau fc Pozzetti|[2000D . 

Finally, we comment on the possibility of the resid- 
ual light from within the Galaxy in the detected COB. 
In Section I2.3| we show that the effect of the residual 
ZL is little, since the systematic difference between the 
two data subsets taken at the systematically different 
heliocentric distances is no larger than those from other 
factors, which are actually include d in our error budget . 
It is consistent with the finding of Manner et al.l ()1974[ ) 
that there are no detectable ZL beyond R = 3.3 AU. 



While the starlight subtraction at the fainter magnitudes 
than the GSC-II detection limits rely on the star-count 
model down to 32.0 mag, absolute contributions of these 
faintest stars are less than ~0.1 bgu at the both bands at 
the Galactic latitudes |6| > 60°. Therefore, their uncer- 
tainties cannot significantly affect our results. Differen- 
tial starlight brightness at the faint limit of 32.0 mag is 
less than 1 x 10 -5 bgu mag^ 1 , hence even fainter stars 
are not likely to make considerable contribution to the 
detected COB levels. The DGL components which do 
not correlate with 100 /im emission are other sources of 
radiation that cannot be removed from the IPP bright- 
ness in the present analysis. Line and continuum emis- 
sions from the interstellar gas might be candidates of 
such components, if their spatial distributions are signif- 
icantly different from that of the interstellar dust pro- 
ducing the dominant component of the DGL. However, 
no such emissions are currently known at the level of ex- 
ceeding a few times 0.1 bgu (see Section |33|) . In any 
case, the residual light from within the Galaxy in the 
detected COB would indicate the true COB brightness 
even smaller than the present estimates, and hence fur- 
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ther strengthen our conclusion that the bulk of the COB 
has already been resolved into galaxies by the current 
deepest observations. 

6. SUMMARY 

We present the results from an analysis of the Pioneer 
10/11 IPP measurements. The used data were obtained 
when the spacecrafts were at the heliocentric distances 
R > 3.3 AU, where brightness of the ZL is below the 
detection limit of the instruments. We carefully examine 
the data quality and select about 70% of the whole data 
as usable for the present analysis. The selected data are 
integrated into sky maps of the high Galactic latitudes 
(\b\ > 35°) at the two bands of the IPP, B IPP (-0.44 ^m) 
and i?ipp (~0.64 /im). Systematic and random uncer- 
tainties of the measurements are estimated to be within 
1% and 2%, respectively, at the both bands. 

Galactic starlight in the IPP brightness is subtracted 
using the two all-sky star catalogs, namely, the Tycho- 
2 catalog and the HST GSC-II. Contribution of fainter 
stars than the detection limits of these catalogs are esti- 
mated from the stellar population synthesis code TRILE- 
GAL, while they are found to be negligible. The DGL is 
further separated out from the IPP brightness by making 
use of its correlation with thermal 100 fim emission taken 



from the IRAS-COBE/DIRBE map. The final resid- 
ual light represents the extragalactic component, i.e., the 
COB. 

The derived COB brightness is 1.8 ± 0.9 and 1.2 ± 
0.9 bgu at £?ipp and i?ipp, respectively, or 7.9 ± 4.0 and 
7.7 ± 5.8 nW m~ 2 sr _1 . From comparison with the in- 
tegrated brightness of galaxies, we conclude that bulk of 
the COB is comprised of normal galaxies which have al- 
ready been resolved in the current deepest observations. 
On the other hand, the derived mean DGL-to-100 /im 
brightness ratios, 2.1 x 10~ 3 and 4.6 x 10~ 3 at the two 
bands, are roughly consistent with the previous observa- 
tions toward the denser dust regions. The red color of 
these ratios confirm the presence of ERE in the diffuse 
ISM. 
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